COMPUTATION OF POLARIZED METRIZED GRAPH INVARIANTS 
BY USING DISCRETE LAPLACIAN MATRIX 



ZUBEYIR CINKIR 

Abstract. Several invariants of polarized metrized graphs and their applications in 

■ Arithmetic Geometry are studied recently. In this paper, we give fast algorithms to 
. compute these invariants by expressing them in terms of the discrete Laplacian matrix 

■ and its pseudo inverse. Algorithms we give can be used for both symbolic and numerical 
. computations. We present various examples to illustrate the implementation of these 

'"^ ■ algorithms. 

(N 

^ . 1. Introduction 

\ Let X be a geometrically connected curve of genus g > 2 over a field k. Suppose 

■ A; is either a number field or the function field of a smooth projective curve Y over a 
field. Assume that X has a semistable model X over SS, where SS = SpecOk if A; is a 
number field and SS = y if A; is a function field. Let N{v) be the local factor related 
to the product formula for k. In this context, the self intersection of the relative 
dualizing sheaf ojx/ss^ is an important quantity both in geometric and arithmetic case. 

Tij- ■ Lower and upper bounds to this quantity are of interest in diophantine geometry. 

^ ! In 1993, Zhang [2l] expressed uj^ in terms of u^, the self intersection of the admissible 

^ I relative dualizing sheaf Ua associated to a minimal regular model of X: 

§■ (1) ^l/ss = ^l + ^4Xv)\ogN{v), 

where v runs over the set of non-archimedean places of Ac, and e{X^) is a certain local 
^ ! invariant of the reduction graph R{X^) associated to the completion of X at a place v of 

k. 

In 2010, Zhang [221 Corollary 1.3.2] expressed a;^ in terms of (A^, A^), the arithmetic 
self intersection (equal to the canonical height) of the Gross-Schoen cycle A^ C X x X x X: 

(2) = |^(^^ + 1^ E ^(^^) ^(-)' 

where v runs over the set of places of k, and if{Xy) is a certain local invariant again linked 
to the reduction graph i?(X„) at each place v of k. 

On the other hand, we have the following equality for a semistable fibration / : X — y 
SS by Nother's formula: 

(3) ^x/ss = 12 deg (ux/ss) -^S{X^), 
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where d{X^) is the total number of the singular points in the fiber over v. Moreover, 
Zhang p5] showed that 

(4) degMu;,/ss) = ^^-^{A^,A^) + Y,KX.)\ogN{v), 

where v runs over the set of places of k, and X{X^) is another local invariant associated 
to the reduction graphs R{X^). We refer [14j for other connections between X{Xy) and 
some invariants of complex moduli space of curves of genus higher than 1. 

Whenever f is a non-archimedean place, the local invariants ip{X^), e{Xy), X{Xy) are 
defined as the invariants of the corresponding polarized metrized graph. Note that the 
reduction graph R{X) of any semistable curve X of genus g over a discrete valuation ring 
is a polarized metrized graph of genus g. 

We also have invariants Z{r), T(r) and 6{r) (see §2] below) of a polarized metrized 
graph r, which are closely related to the local invariants given above. 

Explicit computations of these local invariants were done only for some curves of genus 
less than or equal to 4 (see |I7], [18], [H], [21], [22], [23], [I3], [5]). Other than some 
families of polarized metrized graph of certain types, explicit computations become a 
huge intricate task for a person as soon as g gets larger than 3. Thus, one needs a 
computer algorithm to do such computations. We have such an algorithm to compute 
r(r) (see [5] and [H]). Also, X. Faber gave an algorithm to compute the invariants e(r), 
^(T), X(T) and Z{T) (see [U] and [I2]). 

In this paper, we provide a fast computer algorithm that can be used for both sym- 
bolic and numeric computation of each of the polarized metrized graph invariants. This 
algorithm is faster, because it takes advantage of two important facts: All of the effective 
resistance computations can be handled by just one matrix inversion (see Equation ffTTl) 
and Lemma [3.21) . and the computations of the invariants e(r), v^(r), A(r) and Z(r) can 
be reduced to the computations of t(T) and 9{T) (see Theorem 12.31) . 

In §21 we give a short revision of a polarized metrized graph and definitions of its invari- 
ants (see Lemma \2.1\ Equation (|8]) and Equation ([9])). Using Theorem 12. 3[ the problem 
of computing the invariants e(r), f(T), X(T) and Z(T) are reduced to the computations 
of r(r) and 6'(r). Using these connections and our previous results about t(T) and 0(T), 
we explain how to deal with self loops and multiple edges, if desired, as the initial step of 
the algorithm. 

In ^ we define discrete Laplacian matrices associated to polarized metrized graphs. 
We express 9{r) in terms of the discrete Laplacian matrix L and its pseudo inverse L"*" 
(see Theorem 13. 6p . Then we use our previous result on riT) to express the invariants 
e(r), v^(r), A(r) and Z(r) in terms of the entries of L and L+. This gives us the main 
result in this paper. Namely, each of these invariants can be symbolically or numerically 
computed with the algorithm we provide. 

Finally, we give various examples to show the implementation of our algorithm in ^ 

2. Polarized Metrized graphs and their invariants 

In this section, we first give brief descriptions of a metrized graph F, a polarized metrized 
graph (F,q), invariants t(F), 0(T), e(F), (p(T), A(F) and Z(T) associated to (F,q). 

A metrized graph F is a finite connected graph equipped with a distinguished parametriza- 
tion of each of its edges. A metrized graph F can have multiple edges and self-loops. For 
any given p G F, the number v{p) of directions emanating from p will be called the valence 
of p. By definition, there can be only finitely many p G F with v{p) 7^ 2. 
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For a metrized graph F, we will denote a vertex set for F by ^(F). We require that 
V{r) be finite and non-empty and that p G ^(F) for each p G F if v{p) ^ 2. For a given 
metrized graph F, it is possible to enlarge the vertex set ^(F) by considering additional 
valence 2 points as vertices. 

For a given metrized graph F with vertex set V{r), the set of edges of F is the set of 
closed line segments with end points in V(T). We will denote the set of edges of F by 
E(T). However, if is an edge, by F — Cj we mean the graph obtained by deleting the 
interior of e,. 

We define the genus of F to be the first Betti number g{r) := e — v + 1 of the graph F, 
where e and v are the number of edges and vertices of F, respectively. 

We denote the length of an edge e, G E(T) by Lj, which represents a positive real 
number. The total length of F, which is denoted by -^(F), is given by ^(F) = Yl'i=i 

If a metrized graph F is viewed as a resistive electric circuit with terminals at x and y, 
with the resistance in each edge given by its length, then r(x, y) is the effective resistance 
between x and y when unit current enters at y and exits at x. 

For any x, y in F, the resistance function r{x,y) on F is a symmetric function in x 
and y, and it satisfies r{x,x) = 0. For each vertex set ^(F), r{x,y) is continuous on 
F as a function of two variables and r{x,y) > for all x, y in F. For proofs of these 
facts and connections between the resistance function and the voltage function on F, see 
articles [TU], [21 sec 1.5 and sec 6], and [24, Appendix]. The resistance function r{x,y) on 
a metrized graph were also studied in the articles [I] and [B] . 

We will denote by Ri{T), or by Ri if there is no danger of confusion, the resistance 
between the end points of an edge of a graph F when the interior of the edge is 
deleted from F. 

The tau constant t(F) of a metrized graph F was initially defined by Baker and Rumely 
[2], Section 14]. The following lemma gives a description of the tau constant. In particular, 
it implies that the tau constant is positive. 

Lemma 2.1. [21 Lemma 14.4] For any fixed y in T, r(F) = j J-p (^-^r{x,y)Ydx. 

One can find more detailed information on r(F) in articles [5], [6], [7] and [9]. 

On a metrized graph F, we have a canonical measure /ican fii'st studied by Chinburg and 
Rumely [TD]. See the articles [2] and [6J for several interpretations of fican- The following 
theorem gives an explicit description of fican- 

Theorem 2.2. [lOt Theorem 2.11] For a given metrized graph F, let Li and Ri be defined 
as before. Then we have 



/^can(a;) = ^ {l-^v{p))5p{x) + 



dx 



- Li + Ri 

Let F be a metrized graph and let q : F — )■ N be a function supported on the set of 
vertices of F. The canonical divisor K of (F, q) is defined to be the following divisor on 
F: 

(5) ^= Yl (^(^') - 2 + 2qb))p, and 5k{x) = ^ {v{p) - 2 + 2q(p))5p(a;). 
pev(r) pev(r) 



The pair (F, q) is called a polarized metrized graph (pm-graph in short) if q is non-negative 
and K is an effective divisor. Whenever q = 0, (F,q) is called a simple pm-graph. The 
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genus g(T) of a pm-graph (F, q) is defined to be 

(6) 9{r) = l + ^degK = g{T) + ^ q(p). 

pey(r) 

We will simply use g instead of g{T) when there is no danger of confusion. 

Let fj.ad{x) be the admissible measure associated to K (defined by Zhang [211 Lemma 
3.7]). We have 

^ peV{T) e,ei?(r) * 



Then, if we use Theorem \2.2\ we can relate fiadi^) and ficani^) as follows: 

(7) fJ^adix) = ^{2flcan{x) + Sk{x)). 

Moreover, 5k{T) = deg{K) = 2g - 2, and we have ^caniX) = 1 = fJ-adO^)- 

On a pm-graph (F, q), we defined and studied ([5] and [8]) the invariant 9{T) as follows: 

(8) OiT):= Y (^(p)-2 + 2q(p))(t>(g)-2 + 2q(g))r(p,Q). 

P,9ev(r) 

We have ^(F) > for any pm-graph F, since the canonical divisor K is effective. 

Now, we can give definitions of the invariants e(F), v'(F), A(F) and Z{r) (c.f. [251 
Section 4.1]) of F: 

e(r) = / / r{x,y)6K{x)l^ad{x), Z(F) = i / / r{x,y)fiad{x)Had{y), 
^g^j J JrxF ^ J Jfxt 

V.(F) = 3-g ■ Z{T) - i(6(F) + ^(F)), A(F) = ^|-^V^(r) + ^(e(F) + i?(F)). 

We can express invariants given in Equation ([9]) in terms of t(F) and 6(r) ([8, Propo- 
sitions 4.6, 4.7, 4.9 and Theorem 4.8]): 

Theorem 2.3. Let (F, q) be a pm-graph. Then we have 

(5^-2)r(F) , ejT) £(F) (2^ - l)r(F) , g(F) 

g 4g 4 g^ Sg^ 

^ (3^-3)r(F) g(F) (^ + l)f(F) ^ (Ag - 4)r(F) ^ 

^ ' Ag + 2 16^ + 8 16^ + 8 ' ^ ^ g 2g ' 

Remark 2.4. Given a pm-graph (F, g) anc? s G T — V(T), we have q{s) = andv{s) = 2. 
If we enlarge the vertex set by considering s as a vertex with the same q value on s, 0(T) 
does not change. We recall that this process also does not change the value of t(T) (the 
valence property of t(T), see [6], Remark 2.10]j. Conversely, if q{s) = and v{s) = 2 for 
some s G V^(F) with V{r) has at least two elements, then removing s from the vertex set 
of r does not change t(F) and 6{r) (such vertices are called eliminable vertices in [2T1 
pg. 152]j. Using these observations along with Theorem \2.3l we note that e(F), Z{r), 
<f{T) and A(F) do not change under this process. That is, each of these invariants has the 
valence property. 

Let (F, q) be a polarized metrized graph containing a self loop of length L at a vertex p 
with v{p) > 3. Let /3 be a metrized graph obtained from F by deleting this self loop. We 
still have p G V{f3) but its valence is reduced by 2. Now, we consider a function q : /3 — )■ N 
supported on the set of vertices of (3 that satisfies q(p) = q(p) + 1 and q(x) = q(x) for 
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every x E (3 — {p}. Note that (/?, q) is a polarized metrized graph, and that g{T) = g{f3), 
e(T) = i{/3) + L, eiT) = 9{/3) and t(T) = r(/3) + ^ by the additive property of t(T) P 
pg. 15] and by the fact that a self loop of length L has the tau constant value j^. Using 
these facts along with Theorem I2.3[ we obtain the following equalities: 



(10) _ _ _ 

A(r) = Xm + e(r) = 6(/3) + ^L. 

Successive application of the process above gives a pm-graph 13 that is either a graph with 
one vertex and a self loop at this vertex or a graph without self loops. For the first case, 
we can use the following proposition (part of which is nothing but [251 Proposition 4.4.3]. 
Here we give a new proof.): 

Proposition 2.5. Let (F, q) be a pm-graph with one vertex and e > 1 self loops at its 
vertex. Then 



Proof. If F has just one self loop, then r(F) = and 6{T) = 0. Thus, the result follows 
from Theorem 12.31 in this case. Applying the procedure given in Equation ffTOj) . one can 
show that the result still holds if F has more than one self loop. □ 

Hence, we can focus on pm-graphs with no self loop. 

3. Discrete Laplacian 

In this section, we first describe the discrete Laplacian matrix of a finite weighted graph 
with no self loops and multiple edges. Then following [9], we define discrete Laplacian 
matrices associated to a metrized graph F and a pm-graph (F,q). Then we express 
invariants of a pm-graph in terms of the associated discrete Laplacian matrices, their 
pseudo inverse and values of q. This enables us to achieve our main goal in this paper. 
Namely, we derive a feasible algorithm for both numeric and symbolic computations of 
pm-graph invariants. 

To have a well-defined discrete Laplacian matrix L for a metrized graph F, we first 
choose a vertex set V^(F) for F in such a way that there are no self-loops, and no multiple 
edges connecting any two vertices. This can be done by enlarging the vertex set by 
considering additional valence two points as vertices whenever needed. We call such a 
vertex set V{r) adequate. If distinct vertices p and q are the end points of an edge, we 
call them adjacent vertices. 

Let F be a metrized graph with e edges and an adequate vertex set V^(F) containing v 
vertices. Fix an ordering of the vertices in V^(F). Let {Li, L2, ■ ■ ■ , L^} be a labeling of 
the edge lengths. The matrix A = (apg)^xt) given by 

if p = q, or p and q are not adjacent, 
if p 7^ g, and an edge of length connects p and q. 

is called the adjacency matrix of F. Let D = diag(cipp) be the v x v diagonal matrix given 
by dpp = Xlsevcr) ^ps- Then L := D — A is called the discrete Laplacian matrix of F. That 
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is, L = {lpg)yx^ where 

{0, ii p ^ q, and p and q are not adjacent. 

— a p q, and p and g are connected by an edge of length . 
~ J2sev{r)-{p} hs-> p = q 
We define the discrete Laplacian matrix corresponding to pm-graph (F, q) as the dis- 
crete Laplacian L corresponding to the metrized graph F. Important thing is that if the 
vertex set of (F, q) is needed to be enlarged to make it an adequate vertex set we set q 
value for those new vertices added because of self loops or multiple edges. 

Although L is not invertible, it has generalized inverses. In particular, it has the pseudo 
inverse L"*", also known as the Moore-Penrose generalized inverse. The pseudo inverse L"*" 
is uniquely determined by the following properties: 

i) LL+L = L, in) (LL+)^ = LL+, 

ii) L+LL+ = L+, iv) (L+L)^ = L+L. 

For a discrete Laplacian matrix L of size v x v, the following formula for L"*" (see [201 
ch 10]) is crucial for our computations: 

(11) L+= (L--J)"' + -J. 

where J is of size v x v and has all entries 1. 
L and L"*" are symmetric matrices. 

Remark 3.1. We have EpGy(r) ^pg = = EpGy(r) ^pq' f^r each q G ^(F). 

Lemma 3.2. [3], [1], [151 Theorem A] Suppose T is a graph with the discrete Laplacian 
L and the resistance function r{x,y). For the pseudo inverse we have 

r{p, q) = Ipp - 2/+ + /+ , for any p, q e V{T). 

It is important that t(F) can be expressed in terms of the discrete Laplacian matrix 
and its pseudo inverse: 

Theorem 3.3. [Qj Theorem 4.10] Let L be the discrete Laplacian matrix of size v x v for 
a metrized graph T, and let L+ be its Moore-Penrose pseudo inverse. Suppose Pi and qi 
are end points of Ci G E{V). Then we have 

t(F) = -— V L-„-(— + C -2/+ V L,/+J+ + -tr(L+). 

eiGE(r) P^"?' q,s&V(V) 

Moreover, X]g,sev(r) hs^qq^ts = ~ J2e,eE{r) ^p^q^{K^p^ ~ ^Xq) ■ 

Similarly, one can express ^can in terms of the discrete Laplacian matrix and its pseudo 
inverse: 

Proposition 3.4. [9, Prop. 4.12] For a given metrized graph T, let L be its discrete 
Laplacian, and let L"*" be the corresponding pseudo inverse. Suppose Pi and qi are end 
points of ei G E{r). Then we have 

/ican(2^) = XI ^ ^-^(P)) ^pi^) - Yl + ^^qS^P^P^ " ^^ftg, + CJ)^^" 

p€V{r) e,e-E(r) 

Using Equation ([7]) and Proposition 13. 4[ we express [ladix] in terms of the entries of L 
and L"*": 
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Proposition 3.5. For a given pm-graph (F, q), let L be its discrete Laplacian, and let L"*" 
he the corresponding pseudo inverse. Suppose pi and qi are end points of ei G E{T). Then 
we have 

fiadix) = ^ Yl '^P^^) 'I ^ ^^"''^ ^ ^P^i^^^rHP. - '^Km + ^^gJ)"^^- 

The next important observation is that 0{T) can be expressed in terms of L and L"*": 
Theorem 3.6. Let (F, q) he a pm-graph. Then 

e{V) = 2{2-g-2) {v{p)-2 + 2q{p))l%-2 ^ {v{p) + 2q{p)){v{q) + 2q{q))i;^. 

Proof. Since deg{K) = 2g — 2, ^sgy(r)(''-'('S) ~2 + 2q(s)) = 2g — 2. Then the result follows 
from Equation (|H]), Lemma [3.21 and Remark 13.11 □ 

Whenever q = on F, i.e. F is a simple graph, we have g = g = e — v + 1 and 

(12) ^(F) = 2{2g -2) Y (<P) " '^)Kp " 2 E <p)<<D^p,- 

p(iV{T) p,q&V{V) 

Moreover, if F is both simple and r-regular, we have e = 1^;, g = ^-^v + 1 and we have 

(13) ^(F) = 2i;(r-2)V(L+), 
since ^pgy(r) Ipq = 0. 

Given a pm-graph (F, q), we can give the steps of the algorithm to compute the invari- 
ants e(F), v^(F), A(F) and Z{r) as follows: 

Set L be the discrete Laplacian matrix corresponding to F. (Choose an adequate vertex set so that 
there will be no self loops or no multiple edges. Alternatively, use Equation PU)) for self loops.) 
Set V be the number of rows in L (i.e., the number of vertices in F). 

Set e be the number of nonzero entries above the diagonal of L (i.e., the number of edges in F). 
Set ^(F) be the negative of the sum of the reciprocals of nonzero entries above the diagonal of L (i.e., 
the total length of F). 

Compute g by using Equation 

Compute L+, pseudo inverse of L, by using the formula given in Equation (fTTj) . 
Compute t{T) by using Theorem 13.31 

Define a function v whose value at i is, v{i), the number of nonzero off diagonal entries in the z-th row 
of L (i.e., v{i) is the valence of the vertex i). 
Compute 9(r) by using Theorem 13.61 
Compute e(F), (^(F), A(F) and Z(F) by using Theorem[ 



Note that the most costly part of this algorithm is to obtain L"*" from L. 



4. Examples 

In this section, we give various examples illustrating the implementation of the algo- 
rithm described in the previous section. 

A symbolic computation is exemplified as follows: 
Example 1: 

Let F be a simple pm-graph as illustrated in Figure [H In this case, we have £(F) = 
^ + l + ^ + 2~^'e~^l^ discrete Laplacian matrix corresponding to F is given as 
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Figure 1. Complete graph with vertices {p,q,s,t} and edges of lengths 



111 

' e 



- -I 



follows: 



L 



a + b + c 


—a 


-b 


—c 


—a 


a + d + e 


-d 


—e 


-b 


-d 


b + d + f 


-f 


—c 


— e 


-f 


c + e + f 



We first compute L"*" by using Equation fllip . then use the algorithms given above to 
derive the following results: 

1 A + 2B 



TiT) = -m - 
m = Im - 



z(r) 



6C ' 
2A + 7B 

9C ' 
A + 8B 



9iT) 

m 



6A + 8B 
C ' 

28^^^^ + 28C' 



5A + AB 
9C ' 



108"'"' lose 

where A = ab + ac + be + ad + bd + ae + ce + de + bf + cf + df + ef, B = cd + be + af, 
C = abc+acd+bcd+abe+bce+ade+bde+cde+abf+acf+adf+bdf+cdf+aef+bef+cef. 

The following example has mostly numeric computations and some symbolic computa- 
tions because of non-zero q. 

Example 2: 

Let r be a pm-graph (F, q) such that F is the complete graph on 4 vertices {p, q, s, t} 
as in Figure [H Suppose that each edge of F has length | and that q(p) = q(g) = q(s) = 
q(t) = k for some nonnegative integer k. In this case, we have £(F) = 1, the following 
discrete Laplacian matrix and its pseudo inverse: 



1 
-6 
-6 
-6 



8 -6 



-6 
-6 



-6 
-6 

.8 

-6 1^ 



-6 
-6 
-6 



1 


-1 


-1 


-1 " 


32 


96 


96 


96 


-1 


1 


-1 


-1 


96 


32 


96 


96 


-1 


-1 


1 


-1 


96 


96 


32 


96 


-1 


-1 


-1 


1 


96 


96 


96 


32 . 



Then we obtain these results: 
5 

96' 

16A;2 + 42A; + 25 



r(F) 
A(F) 



^(F) 
Z(F) 



{l + 2kf, 
48A;2 + 88A; 



37 



96(4A; 

17 



3)^ 



96P + loOfc + 17 
96(4A; + 3) ' 
{12k + 11)(2A; + 1) 
12(4A; + 3) ' 
37 ^/■p^ _ 11 



32(8A; + 7) 

In particular, if A; = 0, ^(F) = 1, (^(F) = ^, A(F) = ^, Z(F) = ^ and e(F) = 
Sometimes we are given pm-graphs containing self loops or multiple edges. We exem- 
plify our strategy for such cases in detail as follows: 
Example 3: 
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d+e 




Figure 2. A pm-graph with a self loop, nonzero q and two multiple edges. 



In this example, we consider a pm-graph (r,q) as the graph I in Figure Ej where the 
edge lengths and the values of q are illustrated. In this case, F has a self loop and two 
multiple edges. We first ignore the vertex with q value and valence 2. To avoid having 
multiple edges, we add a vertex to one of the multiple edges such that the added vertex has 
q value 0. We can deal with the self loop by either applying Equation ffTOl) or considering 
two additional vertex with q value on the self loop. These cases are illustrated by graphs 
II and III in Figure O Thus, we can compute polarized metrized graph invariants of F 
in two different ways, and the results are as follows: 

_ t + 2c + 3(a + ^) ^ 4806c 

^ ' 12 ' ^ ' h + 2c 

11a 1162 + 764fec + 44c2 , .-ps _ 9a ?,{h'^ + 2Ahc + Ac^) 21 A 

'^^^^-^+ 72(6 + 2c) ^T' ^^~25^ 25(6 + 2c) ^ 

11a 1162 + 7646c + 44c2 23a 236^ + 8126c + 92c2 3A 

'^^^"l2"+ 36(6 + 2c) ^^^^" 576 + 1728(6 + 2c) + 16 ' 

where A = d + e, £(F) = 3a + 6 + 2c + c? + e, ^(F) = 2 and ^(F) = 12. 

The following example is about the computation of invariants of a class of simple pm- 
graphs. It is given to show both symbolic computations and different aspects of numerical 
computations which would be critical during the implementation of our algorithm. 

Example 4: 

A ladder graph L„(a,6) is planar graph that looks like a ladder with n rungs. It has 
2n vertices and 3t2, — 2 edges. Edges looking like rungs are of lengths 6, and each of 
the remaining 2{n — 1) edges has length a. Thus, £{Ln{a,b)) = 2{n — l)a + nb and 
g{r) = n — 1. Figure [3] shows an example. Using the same notation for the corresponding 
simple pm-graph, we obtain the results given in the tables below: 

We used Mathematica [19] to do the symbolic computations given in Tables [21 [3] and 
m As the number of vertices gets larger the expressions in the results become more com- 
plicated, so we listed the results only for small values of n. Again we used Mathematica 
to do exact computations given in Table O 
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a a 
b 



a a 
b 



Figure 3. Ladder graph L3{a,b). 



n 


c 


V 


g{Ln{a,b)) 




2 


4 


4 


1 


2(a + b) 


3 


7 


6 


2 


4a + 36 


4 


10 


8 


3 


2(3a + 2b) 


5 


13 


10 


4 


8a + 56 



Table 1. Number of edges, vertices, genus and length of Ln{a,b) for n e {2, 3, 4, 5}. 



n 


r(r) 


0(r) 
«r)) 


2 


1 

12 





3 


8(!2^11„6^7/,2 
12(2a-r;-!/j)(la-r3fc) 


26(2«+M 
(2a+36)(la+36) 


4 


12a-i +36a^b+3Sa'^ b'^+lSab^ +36'* 


4a4 +20a^ b+26a2 b"^ +10ab^ +b'^ 


24(a+6) (3a+26) (2a2 +4a6+fe2 ) 


(a+b) (3a+26) (2a2 +4a6+&2 ) 


5 


1 28a-"'' +496a'' 6+704a-^ 6^ ^4yg„2 _|_ gQ^54 ^ ig^s 


2(128(/'+656a*6+1104(r^62+716a2fc» + 160a6'*+965-) 


12(8a+56) (4a2 +6a6+62 ) (4a2 + 10a6+5fe2 ) 


(8a+56) {4a2 +6a6+b2) (4a2 +10a6+562 ) 



Table 2. For pm-graph F = L„(a, 6), r(r) and ^(F) for n e {2, 3, 4, 5}. 



n 


^(r) 


A(r) 
^(r) 


2 





1 

12 


3 


2a2+2a6+62 
3(2a+3fe)(4a+36) 


1 

10 


4 


72a'' +192a'^6+ 164a2 6^ +60a6^ +96"* 


5a+3b 


72{a+6) (3a+26) (2a2 +4afe+62 ) 


14(3a+2b) 


5 


(2a+6) ( 16a'' +52a3 6+51a2 6^ + IGab^ +26*) 


5(2a+6) 


(8a+5b) (4a2 +6ab+b^ ) (4a2 + lOab+Sb^ ) 


9(8a+56) 



Table 3. For pm-graph F = L„(a,6), ip{r) and A(F) for n e {2, 3, 4, 5}. 



n 


^(r) 
^(r) 


z(r) 
^(r) 


2 





1 

12 


3 


4a2+10afe+562 


{a+bf 


3(2a+3fe)(4a+3fe) 


2(2a+3fe)(4a+3fe) 


4 


36a'* +132a3 6+154a2 b'^+66ab^ +96"* 


36a4 + 120a3fe+ 134a2 b'^+60ab^ +9b'^ 


18(a+6) (3a+26) (20^ +4ab+b'^ ) 


108(a+6) (3a+26) (2a2 +4a6+62 ) 


5 


64a5 +288a* 6+452a3 6^ _|_298a2 6^ +80a6* +76^ 


5(2a+b)(16a4 +600^6+ 73a262+32a63+464) 


(8a+5fi) (4a2 +6a6+62 ) (4a2 + I0ab+5b2 ) 


24(8a+56) (4a2 +Qab+b^ ) (4a2 +10a6+562 ) 



Table 4. For pm-graph F = Ln{a,b), e(F) and Z(F) for n e {2, 3, 4, 5}. 
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lb 


r(r) 


m 
m 


v(r) 


A(r) 

^(r) 


m 


Z{T) 

m 


5 


661 


5546 


411 


5 


1189 


925 


10868 


2717 


2717 


39 


2717 


21736 


10 


2107 
37829 


554308 
37829 


30329 
71676 


15 
76 


344578 
340461 


210215 
6128298 


15 


3061011619 
56529128700 


180955287578 
4710760725 


2384321993 
3411240525 


665 
2494 


155613041207 
98925975225 


2950668709 
92330910210 


20 


105284865781 
1971566979888 


6020905705851 
82148624162 


12183994532757 
12486590872624 


380 
1131 


6652614900537 
3121647718156 


10979128575725 
355867839869784 



Table 5. For pm-graph F = L„(l,l), computations of t(F), ^(F), t/?(F), 
A(F), e(F) and Z{T) when n e {5, 10, 15, 20}. 



n 


r(r) 

m 


eir) 

m 


v(r) 

m 


A(r) 
^(r) 


m 


z{r) 

m 


500 


0.05134130 


55155.801168 


27.6397 


7.00233 


55.4713 


0.0278941 


1000 


0.05129978 


221422.325273 


55.4174 


13.9468 


111.027 


0.0278359 


5000 


0.05126659 


5551550.777737 


277.639 


69.5023 


555.471 


0.0277894 


10000 


0.05126245 


22214198.121181 


555.417 


138.947 


1111.03 


0.0277836 


15000 


0.05126107 


49987945.991478 


833.194 


208.391 


1666.58 


0.0277816 


20000 


0.05126038 


88872829.563811 


1110.97 


277.835 


2222.14 


0.0277807 


25000 


0.05125992 


138868444.360781 


1388.75 


347.279 


2777.69 


0.02778 



Table 6. For pm-graph F = L„(l,l), computations of t(F), 6'(F), v^(F), 
A(F), e(F) and Z(F) when n e {500, 1000, 10000, 15000, 20000, 25000}. 



Finally, we used Matlab [16], which does machine arithmetic, to obtain the results in 
Table El Note that when n = 25000 in Table El L„(l, 1) has 50000 vertices and 75000 - 2 
edges. As such computations would be possible on a computer with high memory and 
processing speed, we used Mac Pro with processor 2 x 2.93 GHz 6-core Intel Xeon (24 
hyper-threading in total) and memory 32 GB 1333 MHz DDR3 to obtain these results. 

Symbolic computations are clearly the most costly computations. Among the numerical 
computations, we can have exact arithmetic (arithmetic with numbers having an infinite 
number of significant figures, i.e, with numbers having infinite precision), precision arith- 
metic (which involves the numbers with precision more than 18) and arithmetic with 
machine numbers. 
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